Suiting apparatus and suiting method

ABSTRACT

With respect to input data having the spread of errors or the like, “a suiting method” is to obtain a suiting function, which exist within a region of this spread, and which is highly robust. A suiting region (between an upper bound function M + (x i ) and a lower bound function M − (x i )) is set based on the input data with the spread of errors. An inequality system is established with respect to the suiting region and the suiting function having undetermined parameters. A distance from a boundary of a solution region of the parameter values, corresponding to a boundary of the suiting region, is set, and a suiting function f(x) which has the parameter values maximizing the distance, is found. A broken line exhibits a suiting function f′(x) having low robustness, and a solid line exhibits a most-robust suiting function f(x) obtained from an application of this invention.

TECHNICAL FIELD OF THE INVENTION

The present invention relates to an apparatus and a method which can be utilized for a design of, for example, a digital filter and so on, and which gives robustness to a function that may be used for such a design.

BACKGROUND OF THE INVENTION

The “robustness” is defined as a robustness of a behavior of a system against deviations. Since an actual system involves various kinds of deviations, a guarantee of the robustness is a very important issue in engineering. For example, when a behavior of a feedback loop system is not practically influenced by the introduction of a certain level of external disturbances thereinto, it is said that the system has a “robustness”. As for the deviations, there are not only such external disturbances, but also internal variations of some parameter values caused by characteristic changes due to time passing of parts of the system.

The present inventor has already proposed a suiting method, in which a suiting region is previously set with respect to a physical characteristic, and functions (i.e. suiting functions) contained in the suiting region are found as a family of functions with parameters (cf. Unexamined Japanese Patent Publication (HEI) 7-226656)). According to the previous suiting method, a spread solution region is obtained, and, by using the spread solution region, a parameter design could be performed for a family of functions with a robustness against variations of parameter values.

In the parameter design according to a conventional method of optimizing method, for example, by minimizing a target function, only one optimized set of parameter values is obtained. Due to the obtained optimized parameter values being only one set without any spread, in the parameter design, it is impossible to take a quantitative robustness into consideration with respect to a change of a system behavior which is caused by the introduction of deviations of parameter values thereinto.

Also, in the aforementioned previous suiting method, in order to implement a suiting function, corresponding to the obtained parameter values, to a signal-processing device or the like, it is necessary to select only one set of parameter values from the solution region of the parameter values. Nevertheless, in the aforementioned previous suiting method, a general selection method for promoting a robustness of the suiting function is not especially prepared.

Further, in the aforementioned previous suiting method, there is a case where only one solution is calculated rather than a solution region. In order to obtain a solution uniquely, it is necessary to select how solution is calculated by using additional conditions because a solution of the suiting method ordinarily obtains a solution region. Nevertheless, in the past, the setting of uniformizing-conditions was not performed for promoting the robustness of the suiting function. Thus, generally, the obtained parameter values did not have the desired desired robustness.

Accordingly, in view of the above-mentioned problems, the present invention is directed to an improvement of the above-mentioned suiting method, such that solutions with robustness, and therefore robustly-stable physical systems, can be obtained when a solution region or a unique solution is needed.

DISCLOSURE OF THE INVENTION

In spite of a deviation being introduced into parameter values corresponding to the suiting function, when the parameter values still remain in a suiting region, this deviation is defined as an allowable deviation of this parameter values. In a suiting method according to the present invention, a robust parameter value stands for a parameter value which has large allowable deviations, and a most robust parameter value stands for a parameter value which has the largest allowable deviation.

A suiting apparatus according to the present invention at least comprises means for setting a suiting region corresponding to an allowable range of a physical characteristic, and a suiting function determination means for finding at least one suiting function at least approximately contained in a range of the suiting region, by solving simultaneous inequalities. The suiting apparatus is characterized in that the suiting function determination means sets the suiting function as a family of functions with certain parameters, sets a distance from a boundary of a solution region corresponding to a boundary of the suiting region, and finds a suiting function with robust parameter values by solving the simultaneous inequalities under a condition that the distance becomes larger than a beforehand given level or is maximized.

Also, a suiting method according to the present invention at least comprises a first step of setting a suiting region corresponding to an allowable range of a physical characteristic, and a second step of finding at least one suiting function at least approximately contained in a range of the suiting region. The suiting method is characterized in that the second step sets the suiting function as a family of functions with certain parameters, sets a distance from a boundary of the solution region corresponding to a boundary of the suiting region, and finds a suiting function with robust parameter values by solving the simultaneous inequalities under a condition that the distance becomes larger than a beforehand given level or is maximized.

In this specification, “a physical characteristic” stands for a desired characteristic, desired form, measurement value, communication signal, pattern signal, etc. Also, “a suiting” stands for a method of deriving a function (hereinafter referred to as a suiting function) at least approximately within a region corresponding to a spread of an input (hereinafter referred to as a suiting region), and more concretely a method of deriving a set of parameter values of the suiting function.

This suiting method involves an enormous amount of operations even in a case where the suiting function is expressed in a simple form. Accordingly, it is impossible to find the solution by manual calculation, and it only becomes possible by using a digital computer etc. Also, “a suiting method” stands for a method of executing the suiting, and “suiting apparatus” stands for an apparatus for executing the suiting method. Namely, the suiting apparatus stands for a digital computer or the like by which the suiting method is executed.

A typical suiting method will be explained below using a vector notation.

A representative example of a suiting method is expressed as the procedure of setting an appropriate function with parameters as f(x) and finding parameter values satisfying:

M⁻(x)<f(x)<M⁺(x)(∀×εD)  (1)

Herein: D indicates a domain in a multidimensional space. The inequality system (1) is called a determinating inequality system for the suiting, and the f(x) satisfying inequality system (1) is called a suiting function. Also, M⁺(x), which determines the upper bound of the suiting function f(x), is called an upper bound function, and M⁻(x), which determines the lower bound of the suiting function f(x), is called a lower bound function. The entire range determined by upper bound function M⁺(x) and lower bound function M⁻(x) is called a suiting region T(={(x, y)|×εD, M⁻(x)<y<M⁺(x)}). Upper bound function M⁺(x) and lower bound function M⁻(x) may be set as values of +∞ and −∞, respectively. Also, a suiting region may be set by a center value and an allowable width with respect to the center value.

The inequality system (1) is expressed by a finite number of simultaneous inequalities. A range of the parameter values of f(x) which satisfy all of the finite numbers of inequalities is the solution region. This solution region is obtained by transforming the suiting region to a parameter space. Accordingly, one point in the solution region corresponds to a suiting function which attains a value within the suiting region.

When a suiting function is selected from the solution region, a distance from a boundary of a parameter solution region corresponding to a boundary of the suiting region is set, and a suiting function with robust parameter values is selected such that the distance becomes larger than a beforehand given level or maximum.

Note that in the suiting apparatus and the suiting method according to the present invention, the suiting function f(x) is not only a real-value function, but also may be a functional, an operator, etc., and further may be a recursion relations of a sequence. Also, the domain D may be a discrete set of points or a continuous range, which may be disconnected.

Moreover, it is possible that a plurality of suiting functions to be suited is set, if their parameters are common. In this case, respective different suiting regions may be set with respect to functions to be suited.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a view showing an upper bound function, a lower bound function, and a suiting function of Embodiment 1 of the present invention.

FIG. 2 is a view showing a setting of a distance from a boundary of a solution region to a certain parameter value.

FIG. 3 is a view showing a comparison of a maximal robust solution with a solution obtained by a conventional method.

FIG. 4 is a view showing Embodiment 2 of the present invention.

FIG. 5 is a block diagram of the structure of Embodiment 3 of the present invention.

FIG. 6 is a view of an output screen displaying a suiting function obtained by Embodiment 3 of the present invention.

FIG. 7 is a flow chart showing an operation of Embodiment 4 of the present invention.

FIG. 8 is a block diagram of an identification system of Embodiment 6 of the present invention.

DESCRIPTION THE PREFERRED EMBODIMENT

The embodiments of the present invention will now be described below.

First, a suiting apparatus and a suiting method according to the present invention will be explained by Embodiment 1 and Embodiment 2, and then an explanation will be made of the application of the suiting apparatus and method to various systems in the technological and industrial fields by Embodiment 3 and subsequent embodiments.

EMBODIMENT 1

Embodiment 1 is described below.

Here, solution methods for a robust solution, namely robust suiting function, by a suiting method are explained, using the following steps:

1. Reduction of problem of suiting method to inequality system;

2. Solution method of robust solution.

1. Reduction of Problem of Suiting Method to Inequality System

A suiting method will be explained using FIG. 1. In this embodiment, the sequences of points P_(j)=(x_(j), y_(j)) (j=1, 2, . . . , n) are produced based on, for example, input data. The domain D on the axis x, which limits the suiting function becomes {x₁, x₂, . . . , x_(n)}. Considering the spread due to error of the points P_(i), etc. on x=x_(i), upper bound function M⁺(x_(i)) is determined by the upper bound value that can be given by this point P_(i), and lower bound function M⁻(x_(i)) is determined by the lower bound value. Such a set up is carried out with respect to the sequences of points P_(i). Namely, the sequences of points P⁺ ₁=(x₁, y⁺ ₁), P⁺ ₂=(x₂, y⁺ ₂), . . . , and P⁺ _(n)=(x_(n), y⁺ _(n)) corresponding to the upper bound functions, and the sequences of points P⁻=(x₁, y⁻ ₁), P⁻ ₂=(x₂, y⁻ ₂), . . . , P⁻ _(n)=(x_(n), y⁻ _(n)) corresponding to the lower bound functions are determined. Note, y⁺ _(i) is M⁺(x_(i)), and y⁻ _(i) is M⁻(x_(i)). The suiting region T is determined by all of the points of the sequences.

Assume that the suiting function can be expressed by a rational function of the first degree in x as follows:

f(x)=(a₀+a₁x)/(b₀+b₁x)  (∀×εD)  (2)

Herein: a₀, a₁, b₀, and b₁ are real parameters

When the rational function (2) is substituted in inequality system (1) and, at the same time, the denominator is canceled under the condition that the denominator (b₀+b₁x_(i)) is positive, inequality system (1) is transformed into the following equivalent inequality systems (3), (4) and (5):

−a₀−a₁x_(i)+b₀M⁺(x_(i))+b₁x_(i)M⁺(x_(i))>0  (3)

a₀+a₁x_(i)−b₀M⁻(x_(i))−b₀x_(i)M⁻(x_(i))>0  (4)

b₀+b₁x_(i)>0  (5)

So as to express inequality systems (3), (4) and (5) as inner products of vectors, a parameter space of 4-dimensions, corresponding to the number of the parameters, is set. The vectors X, η⁺(x_(i)), η⁻(x_(i)) and η₀(x_(i)) in the parameter space are set as follows:

X=(a₀, a₁, b₀, b₁)  (6)

η⁺(x_(i))=(−1, −x_(i), M⁺(x_(i)), x_(i)M⁺(x_(i)))  (7)

η⁻(x_(i))=(1, x_(i), −M⁻(x_(i)), −x_(i)M⁻(x_(i)))  (8)

η₀(x_(i))=(0, 0, 1, x_(i))  (9)

Herein: X is a parameter vector, η⁺(x_(i)) (hereinafter referred to as an upper bound vector) is a vector representing the limitation of the upper bound function M⁺(x_(i)), η⁻(x_(i)) (hereinafter referred to as a lower bound vector) is a vector representing the limitation of lower bound function M(x_(i)), and η₀(x_(i)) is a vector representing the condition necessary for preventing the rational function (2) from diverging.

When using these vectors, the determinating inequality system of the suiting method expressed by inequality systems (3), (4) and (5) becomes the following inequality systems (10), (11) and (12) expressed by the inner product with the parameter vector X, respectively:

(X, η⁺(x_(i)))>0  (10)

(X, η⁻(x_(i)))>0  (11)

(X, η₀(x_(i)))>0  (12)

In this way, since inequality systems (10), (11) and (12) are given to the input point sequence P_(i) on the domain D, respectively, the problem of finding at least one suiting function is expressed by a finite number of inequalities, and the solution can be actually found. For example, where the number n of the sequences of points P_(i) is 500, there are 500 sets of the inequalities (10), (11) and (12), respectively.

When the finite number of inequalities obtained with respect to all input sequences of points P_(i) are simultaneously exhibited and solved for the parameter vector X, a region, where X exists, is called a solution region S. The solution region S generally becomes a convex cone and is expressed by:

S={s₁X₁+s₂X₂+ . . . +s_(m)X_(m)|s₁, s₂, . . . , s_(m)>0}  (13)

The vectors of X₁=(a₁₀, a₁₁, b₁₀, b₁₁), X₂=(a₂₀, a₂₁, b₂₀, b₂₁), . . . , X_(m)=(a_(m0), a_(m1), b_(m0), b_(m1)) express the vertices of the solution region S. A suiting function f(x), determined by a point X_(s).=(a_(s0), a_(s1), b_(s0), b_(s1)) inside the solution region S, always gives values within the suiting region T.

2. Solution Method of Robust Solution.

Referring to FIGS. 2 and 3, one example of determining an ex_(i)stence of the solution region S and finding the most robust solution X_(s), based on a linear programing method, will be explained. For the sake of simplifying the description, all the vectors of η⁼ (x_(i)), η⁻ (x_(i)) and η₀ (x_(i)) used in the above-mentioned inequality systems (10), (11) and (12) are correspondingly expressed by new vectors Tlil (j=1, 2, . . . , 3n). Thus, inequality systems (10), (11) and (12) are represented by the following equivalent inequality system (14):

(X, η_(j))>0(j=1, 2, . . . , 3n)  (14)

Herein:

η_(i)=η⁺(x_(i))(i=1, 2, . . . ,n)  (15)

η_(i+n)=η⁻(x_(i))  (16)

η_(i+2n)=η₀(x_(i))  (17)

Also, a representation of η_(j) by components is put as (η_(j0), η_(j1), η_(j2), η_(j3)).

In inequality system (14), even if the parameter vector X is multiplied by α, the following inequality system (18) stands:

 (αX, η_(j))=α(X, η_(j))>0  (18)

Therefore, even if X is multilplied by (1/a₀) and the following is set;

X″=(1/a₀)X=(1,a₁/a₀, b₁/a₀, b₁/a₀),

the generality cannot be lost.

Accordingly, the above inequality system (14), expressed by the inner products of the 4-dimensional vectors, is reduced to the following inequality system (19), expressed by the inner products of the 3-dimensional vectors:

0<(X″, η_(j))=η_(j0)+(X′, η′_(j))  (19)

Herein: X′ and η′_(j) are set as follows:

X′=(X′⁽¹⁾, X′⁽²⁾, X′⁽³⁾)=(a₁/a₀, b₀/a₀, b₁/a₀)

η′_(j)=(η_(j1), η_(j2), η_(j3))

A solution region S′ of inequality system (19) is expressed as follows:

S′={s₁X′₁+s₂X′₂+ . . . +s_(m)X′_(m)|s₁, s₂, . . . , s_(m)>0, s₁+s₂+ . . . +s_(m)=1}  (20)

herein: X′_(k) (k=1, 2, . . . , m) is a 3-dimensional vector expressing a vertex of the solution region S′.

A point of the most robust solutions is obtained from the solution region S′. To this end, first, a distance from a boundary of the solution region S′ is set. If j is fixed, inequality (19) geometrically expresses a half-space S_(j) with a boundary which is the following plane H_(j) in R³ (orthogonal coordinates of real 3-components):

 H_(j):ƒ_(j0)+(X′, ƒ′_(j))=0  (21)

Herein: ƒ′_(j) is a normal-direction vector of H_(j); and S_(j) exists on the side of the direction of ƒ′_(j). Since a region, on which inequality (19) stands for any j, is a solution region in R³-space, the solution region S′ is the intersection of each S_(j) (S′=S₁∩S₂∩ . . . ), and the boundary of S′ is the parts of H₁, H₂, . . .

Again, j is fixed. If a vector X′ is in S_(j), and if the distance from H_(j) to X′ is y_(j), the following equation (22) stands:

y_(j)|η′_(j)|=η_(j0)+(X′,η′_(j))(|η′_(j)|=(η_(j1) ²+η_(j2) ²+η_(j3) ²)^(½))  (22)

For the sake of rearrangement, the equation (22) is divided by |η′_(j)| as a whole. Then, if settings of ξ_(j0)=η_(j0)/|η_(j)| and ξ′_(j)=η′_(j)/|η′_(j)| are given, the equation (22) can be expanded as follows:

y_(j)=ξ_(j0)+(X′, ξ′_(j))  (23)

This aspect is shown in FIG. 2 (note, in this drawing, the case of ξ′_(j3)=0 is illustrated for the sake of simplicity). As shown in FIG. 2, the right side region of H_(j) is S_(j), and X′ is on a dashed line parallel to H_(j).

Although the above-mentioned formulation with respect to the distance is established for one fixed j, a distance y from the boundary of S′ to the vector X′ with respect to all of j can be defined by the following:

 y=min (Y₁, Y₂, . . . )  (24)

The formulation implies that a distance from any one of H_(j) to X′ is at least y. In other words, if a distance from X′ to any one of H_(j) is equal to or more than y≦y_(j) stands, and from (23) the following stands:

y≦ξ_(j0)=(X′, ξ′_(j))(j=1, 2, . . . )  (25)

The most robust solution is the farthest solution from the boundary of the solution region S′, because the most robust solution is defined as a set of parameter values with the largest allowable deviations. In inequality system (25), y is a lower limit of the distance from the boundary, thus X′, maximizing y, is farthest from the boundary. Thus, the solution method of the most robust solution is reduced into a linear programming maximizing method by which the subject value y is maximized under the constraint condition of inequality system (25).

The most robust solution X′_(s) and the maximal distance y_(s) represents a disc part of the solution region having a center X′_(s) and a radius Y_(s) Thus, y_(s) implies a quantitative robustness (y_(s) is referred to as a robust radius hereinafter).

In a conventional method, for example, instead of inequality system (25), a linear programming method was used, having the following constraint condition:

y≦η_(j0)+(X′, η′_(j))(j=1, 2, . . . )

Namely, X′, maximizing y under a condition of y>0, was found. Nevertheless, this y does not represent a distance from the boundary, and is merely introduced as an artificial-subject-variable for the purpose of obtaining a solution. Therefore, generally, it is impossible to obtain a robust solution.

FIG. 3 shows a comparison of the conventional solution with the solution obtained by this embodiment. In FIG. 3, S′ is 2-dimensionally illustrated for the sake of simplicity. The solution X′₀ is obtained from the conventional method, and may be out of the solution region S′ even due to a small variation of the parameter values. On the contrary, X′_(s) is the most robust solution, and, even when a deviation is introduced into the parameter values, the set of parameter values still remains in the solution region S′ if a level of the deviation is less than y_(s).

In the above-mentioned embodiment, the unique most robust solution is found, one solution or more than one solution, which have a robustness larger than a predetermined level, may be found (if existing). For example, by setting a relative level of robustness “a” (0<a<1) and by putting a disc neighborhood of the most robust solution X′_(s), with radius y_(s)×a, a solution region with a level of robustness y_(s)×(1−a) may be obtained.

Further, instead of using the maximizing condition for y, by setting an absolute level of robustness y₀, another suiting problem may be established; it makes solve simultaneously inequality system (25) and the following inequality condition that the robustness is more than y₀ as follows:

y₀≦y.

In this case, it is possible to obtain a solution region having a higher level of robustness than y₀. Note, it is possible to suitably determine the above values “a” and “y₀” in accordance with levels of variations of parameter values with respect to an objective system.

From the most robust solution X′_(s)=(X′_(s1), X′_(s2), X′_(s3)) obtained in R³, the following parameter vector is produced:

X_(s)=(α, αX′_(s1), αX′_(s2), αX′_(s3))(α>0)

Thus, by substituting the parameter vector X_(s) in the equation (2), it is possible to obtain most robust suiting functions.

An described above, according to this embodiment, the following effects are obtained:

[1] It is possible to find the most robust suiting functions against variations of the parameter values;

[2] It in possible to quantitatively determine a level of the robustness.

Note, if there is no solution region S′ for the linear programming problem on inequality system (25) or if y_(s)<0, there is no solution region S for the original suiting problem. Namely, by solving a feasible problem of the linear programming problem of inequality system (25), it is possible to determine whether or not the suiting function exists.

In the above-mentioned ezbodant, the formulation (24) of the distance is merely shown by way of example. For example, by using a symmetric positive-definite matrix (metric matrix) gkl(k, l=1, 2, 3), in the parameter space, a “length |Y|” of a vector Y may be defined as follows:

|Y|=(Σ³ _(k,=1)g_(kl)Y_(k)Y_(l))^(½)

and, instead of the equation (22), the distance from the boundary H_(j) to X′ is the following:

y_(j)=(ξ_(j0)+(X′, ξ′_(j)))/|ξ′_(j)|  (26)

Similar to the relationship (25), if X′ maximizing y is found under the constraint condition:

y≦(ξ_(j0)+(X′, ξ′_(j)))/|ξ′_(j)|(j=1, 2, . . . )  (27)

the most robust solution can be obtained in the sense of the weighted distance defined by g_(kl). The definition of the above-mentioned distance is effective in the cases where there are differences among the variations of the parameter values, for example, where the parameter value X₁ varies more greatly than the parameter value X₂.

It is easy to expand the procedures of “1. Reduction of problem of suiting method to inequality system” and “2. Solution method of robust solution” described above with respect to a case where the rational function (2) is generalized to the ratio of linear combinations of two sets of linear independent functions (p₀(x), p₁(x), . . . , p_(M)(x)) and (q₀(x), q₁(x), . . . q_(n)(x)) an the domain D in the n-dimensional space as follows:

f(x)=(a₀p₀(x)+a₁p₁(x)+ . . . +a_(M)p_(M)(x))/(b₀q₀(x)+b₁q₁(x)+ . . . +b_(N)q_(N)(x))

Herein: a₀, a₁, . . . , a_(M), b₀, b₁, . . . and b_(N) are real parameters; x is an n-dimensional vector; and p₀(x), p₁(x), . . . , p_(M)(x), q₀(x), q₁(x), . . . and q_(N)(x) are piecewisely continuous functions on the domain D, such as x^(n), cos x, e^(x) or stop functions, etc.

In the above-mentioned embodiment, with respect to the suiting problem for the ratio of two linear expressions in the parameters, the most robust solution has been found. Even for a general suiting problem being non-linear in the parameters, if a distance from the boundary can be adequately defined, the problem can be reduced into a non-linear optimizing problem for this distance. Also in this case, by solving the optimizing problem, the most robust solution can be obtained.

EMBODIMENT 2

Embodiment 2 of the present invention will now be described. In this embodiment, there is shown a method for selecting the most robust solution from a solution region of parameter values or a part thereof when the region is previously expressed by a formulation similar to the set (20). For the sake of simplicity, S′ of the set (20) is used as the original solution region here. Note, a similar procedure is possible even in a higher dimensional case.

Since S′ exists in 3-dimensional space, in general, its boundary can consist of connected triangles Δ_(j), as shown in FIG. 4. From vertex vectors X′_(a(j)), (j), X′_(b(j)), X′_(c(j)) of each Δ_(j), a vector η′_(j) (≠0), orthogonally intersecting Δ_(j) is calculated. Namely, the vector η_(j), which is directed to the solution region S′, is calculated so as to satisfy the following;

(η′_(j), X′_(b(j))−X′_(a(j)))=(η′_(j), X′_(c(j))−X′_(a(j)))=0.

If y_(j) is set as a distance from X′ to a plane H_(j) containing Δ_(j), the following relationship holds:

y_(j)|η′_(j)|=−(η′_(j), X′_(a(j)))+(X′,η′_(j))  (28)

By letting η′_(j)=−(η′_(j),X′_(a(j))), the equation (28) becomes equivalent to the equation (22). Accordingly, the most robust solution can be obtained in the same manner as in the item “2. Solution method of robust solution” of Embodiment 1.

As is apparent from the foregoing, according to the present Embodiment 2, in the case where the solution region is previously given, the most robust solution can be selected from the solution region of parameter values. Also, according to this method, it is possible to obtain easily the each most robust solution from a same solution region, in accordance with a change of robustness (for example, by changing a distance function).

EMBODIMENT 3

FIG. 5 shows, as Embodiment 3, a schematic structure of an apparatus which executes the suiting method to obtain the most robust solution. The apparatus 11 is, for example, a digital computer, and comprises a memory device 12, a processor 13, an input device 14, an output device 15, and a control device 16, and these elements are mutually connected to each other through bus lines 17 and 18. The memory device 12 includes: a storage area for storing upper bound function Me^(+(x) _(i)), i.e. an upper bound function memory unit 21; a storage area for storing lower bound function M⁻(x_(i)), i.e. a lower bound function memory unit 22; a storage area for storing a program to set a distance from a boundary of a solution region, i.e. a boundary distance setting program unit 23; a storage area for storing a program to obtain a solution maximizing a boundary distance, i.e. a boundary distance maximizing program unit 24; and a storage area for storing a control program, such as an OS, i.e. a control program unit 25. The processor 13 is a CPU or such like. The input device 14 is, for example, a keyboard, a mouse, a numerical file, a digitizer, or a light pen. The output device 15 is, for example, a display, a numerical file, a plotter, or a printer. The control device 16 controls the respective devices for executing the program.

FIG. 6 shows an example of the data displayed on the screen of a display of the output device 15. Upper bound function M⁺(x_(i)) and lower bound function M⁻(x_(i)) are set based on the sequence of points given via the input device 14. A suiting function f(x) found by the suiting method is the following rational function (29) with the denominator and the numerator respectively expressed by fourth-degree polynomials:

f(x)=(A₀+A₁x+A₂x²+A₃x³+A₄x⁴)/(B₀+B₁x+B₂x²+B₃x³+B₄x⁴)  (29)

Herein: A₀, A₂, A₃, A₄, B₀, B₁, B₂, B₃, and B₄ are real parameters.

In FIG. 6, a broken line exhibits a suiting function f′(x) with a small robustness, and a solid line exhibits a most robust function f(x). As is apparent from this, a robustness of the suiting function gives a smoothness to the suiting function per se.

Note, the function (29) represents a squared amplitude characteristic of a digital band pass filter. This will be stated in detail in Embodiment 5.

As described above, according to the suiting apparatus of the present embodiment, the effects similar to those of Embodiment 1 are obtained. Also, the suiting apparatus of this embodiment has a-structure for visually displaying the upper boundary and lower boundary of the suiting region on a screen of the display and, at the same time, displaying the most robust suiting function f(x) on this displayed suiting region in an overlapping manner. Accordingly, a correlation between the suiting function and the spread of the input data can be directly and visually understood, and a past labor to select a suiting function can be saved. Therefore, in a case where the suiting method is applied to, for example, a system design, even if a user does not have detailed knowledge of the object system, it is possible for a user to obtain efficiently a design result, which guaranteed engineering stability, merely by visually estimating the design result and visually instructing changes of specifications and so on.

EMBODIMENT 4

In this embodiment, a procedure for finding at least one suiting function having a minimum degree will now be described.

In general, regarding the same suiting region, a higher-degree suiting function has higher robustness. Thus, it is possible to obtain the minimum-degree suiting function with a desired robustness, by determining an existence of a suiting function, which satisfies the desired robustness, by changing a degree of the function. The robustness is represented by, for example, a robust radius, and is set by a suitable input means (such as the input device 14 in FIG. 5).

One example of the procedure is shown in the flow chart of FIG. 7, and will be described below.

At step 1, an upper bound function and a lower bound function are inputted. At step 2, an initial value of the degree N is made “1”. At step 3, a maximizing problem for the robustness is set and is solved by using a distance function, similar to inequality system (25), which represents a distance from the boundary. At step 4, if a solution, which satisfies the constraint condition, does not exist, the routine proceeds to step 7, and if the solution exists, the routine proceeds to step 5. At step 5, if the robust radius of the obtained solution is smaller than a previously-desired robust radius rot the routine proceeds to step 7, and if larger, the routine proceeds to stop 6. At step 6, a set of parameter values, satisfying the desired robustness, and a robust radius thereof are outputted, and the most robust function is displayed as a graph on the display. Thus, the routine ends.

At step 7, the degree N of the suiting function is incremented by one, because the solution does not exist or does not satisfy the desired robustness. At step 8, if the degree N of the suiting function exceeds the maximum degree Nmax, which can be realized by the apparatus, the routine proceeds to step 9, and, if else, the routine returns to step 3, in which and a most robust function is again searched for. At step 9, a message, for example, “THERE IS NO SOLUTION, SATISFYING DESIRED ROBUSTNESS, OF DEGREE EQUAL TO OR LESS THAN Nmax” is displayed, and thus the routine ends.

Note that, the above-mentioned procedure merely shown a simple algorithm, by way of example, and it is possible to carry out the procedure at a high speed by introducing an algorithm, such as a dichotomy or the like.

Incidentally, the present invention is applicable to various fields, where the suiting method can be used, such as engineering systems represented by transfer functions of analog systems and digital systems in the frequency domain, and engineering systems represented by recursion relations and differential equations allowing non-linear elements and time-variant elements in the time domain, and so on. Also, the present invention is applicable to not only a system design but also an identification of systems or the like.

In all applicable fields and aims, according to the present invention, a most robust suiting function can be similarly obtained, and gives an individually-significant merit to each of the fields and aims. So, two typical applications of the present invention: a design method for an amplitude characteristic of a digital filter in the frequency domain (Embodiment 5), and an identification of systems in the time domain (Embodiment 6) will be described hereinafter.

The design method and the identification method, shown in the following two embodiments, are easily applicable to various systems in engineering, (such as other physical systems including a linear-phase filter, an all-pass filter, analog systems, etc., and signal-processing systems including for example a data-compressing system), only if ways to provide suiting regions are modified using the known techniques. The ways to provide suiting regions of “a linear-phase filter”, of “an all-pass filter”, of “analog systems”, and of “a data-compressing system” are disclosed in the specification of Unexamined Japanese Patent Publication (HEI) 7-226656. Namely, the way to provide a suiting region of “a linear-phase filter” is disclosed in paragraph [0093]; the way to provide a suiting region of “an all-pass filter” is disclosed in paragraphs [0094] to [0098]; the way to provide a suiting region of “analog systems” is disclosed in paragraphs [0099] to [0109]; and the way to provide a suiting region of “a data-compressing system” is disclosed in paragraphs [0110] to 0125], respectively.

EMBODIMENT 5

In this embodiment, a transfer function of a digital filter is designed from a given frequency-amplitude characteristic. For the sake of simplicity, the sampling period T is normalized as 1. In the design, only the amplitude characteristic is specified neglecting the phase characteristic. The domain, on which the characteristic should be specified, is set as D₀:

D₀={0≦Ω≦π}

Although a degree of the transfer function can be arbitrarily set, an explanation will be made of a case where each of the numerator and the denominator of the transfer function is of the 4th-degree.

Transfer function H(z) is defined as follows:

H(z)=(a₀+a₁z⁻¹+a₂z⁻²+a₃z⁻³+a₄z⁻⁴)/(b₀+b₁z⁻¹+b₂z⁻²+b₃z⁻³+b₄z⁻⁴)=P(z)/Q(z)  (30)

The numerator and the denominator are set as P(z) and Q(z), respectively, as described above. The definition (30) relates to the transfer function of an IIR filter, but, in the case of an FIR filter, it is sufficient if Q(z)=b₀.

For the sake of simplicity of calculation, the amplitude characteristic is indicated by a square amplitude characteristic. When substituting z=e^(jω) in H(z), the square amplitude characteristic of the filter is defined by H(e^(jω)) H*(e^(jω)) as follows (* indicates the complex conjugate):

H(e^(jω))H*(e^(jω))=H(e^(jω))H(e^(−jω))=(A₀+A₁ cos ω+A₂ cos²ω+A₃ cos³ω+A₄ cos⁴ω)/(B₀+B₁ cos ω+B₂ cos²ω+B₃ cos³ω+B₄ cos⁴ω)

This square amplitude characteristic is expressed as a function in cos ω as shown above. The numerator and denominator are set as p(cos ω) and q(cos ω) respectively, and f(cos ω) is set as p(cos ω)/q(cos ω). A₀, A₁, A₂, A₃ and A₄ are the quadratic polynomials of a₀, a₁, a₂, a₃ and a₄, and B₀, B₁, B₂, B₃ and B₄ are the quadratic polynomials of b₀, b₁, b₂, b₃ and b₄. Note, the explanation made here uses cos^(k)ω (k=0, 1, . . . ) as a set of the base functions, but the set of the base functions is not necessarily restricted to these. For example, cos kω (k=0, 1, . . . , 4) may be used as another set of the base functions for the suiting on the domain D₀. Upper bound function M⁺ (cos ω) and lower bound function M⁻(cos ω), expressing a desired range of the square amplitude characteristic, are given on the domain D₀, Since the range of cos_ on D₀ is from −1 to +1, letting x=cos ω and D={−1≦x≦1}, the suiting is executed by the following determinating inequality system:

M⁺(x)<f(x)=p(x)/q(x)<M⁻(x)(xεD)

0<q(x)(xεD)

Herein: the positive-valued condition of q(x) is necessary for the stability of the filter. From this, by using the above-mentioned suiting method to obtain the most robust solution, a set of the parameter values A₀, A₁, A₂, A₃ and A₄, and B₀, B₁, B₂, B₃ and B₄ of f(x) are obtained. By using suitable factorizations of the thus obtained square amplitude characteristic p(x)/q(x), the corresponding transfer function P(z)/Q(z) can be obtained.

Note, in the execution of the suiting, generating convex hulls (cf. paragraphs [0048] to [0078] of Unexamined Japanese Patent Publication (HEI) 7-226656) may be utilized to obtain a suiting region with higher precision or more robustness. This embodiment also may be used for an identifuncation of an object system by setting a measured frequency characteristic of this system as a suiting region.

As described above, even if a variation, which has a level corresponding to the robust radius, is introduced into the set of parameter values obtained in accordance with this embodiment, the frequency characteristic of the filter still remains within the suiting region. This implies that, with respect to a design of digital filters, if a level of quantumization noises of the coeffuncients and the input signal is less than the robust radius, a frequency characteristic falls in an allowable range of the specifuncation. Thus, in a case of a design of a fixed-point digital filter, for example, by using the procedure in Embodiment 4 under a condition that a robust radius is more than a level of the quantumization noises, it is possible to design a transfer function, which is of the minimum degree and which is free from deterioration resulting from the quantumization of coeffuncients and input signals.

Also, a similar method can be applied to a design of an analog system under a condition that the robust radius is greater than a level of an aging-deviation or a temperature-deviation of characteristics of each of the parts, whereby the characteristics can be prevented from being subjected to deterioration derived from the deviations.

EMBODIMENT 6

In this eabodiment, a method of system-identifuncation in a time domain will be described. As shown in FIG. 8, in a case where an object system 102 receives a signal x from a signal source 101, and then responds to the signal x with a signal y, in order to identify the object system 102, a model is estimated from the signals x and y by an identifuncation apparatus 103. The signals x and y are discretized, with n-th signals being represented by x_(n) and y_(n), respectively.

When the object system 102 is set, with using suitable functions f_(i), as follows:

S_(y)=Σ^(L) _(i=0)a_(i)f_(i)(y′_(n), x′_(n))

parameters a_(i) are determined, assuming that it is approximately represented by S_(y)=0 Herein: Σ expresses a sum of 0 to L for i, and respective y′_(n) and x′_(n) are set as:

y′_(n)=(y_(n), y_(n−1), . . . , y_(n−M))

x′_(n)=(x_(n), x_(n−1), . . . , x_(n−M))

f₀, f₁, . . . and f_(L) are mutual linearly-independent functions. f_(i)(y′_(n), x′_(n)) may be non-linear in the arguments y′_(n) and x′_(n), and may depend on n (namely, the system is time-variable).

Further, the following definitions are made:

X′=(a₀, a₁, . . . , a_(L)),

η_(n)=(f₀(y′_(n), x′_(n)), f₁(y′_(n), x′_(n)), . . . , f_(L)(y′_(n), x′_(n)))

The respective lower bound function and upper bound function of allowable errors for S_(y) are expressed by M⁻(n) and M⁺(n), which, for example, are derived from measurement errors. In this case, the determinating inequality system for the suiting is as follows:

M⁻(n)<S_(y)(X′, Y′_(n), x′_(n))<M⁺(n)(n=0, 1, . . . )

Then, this is expressed as follows:

M⁻(n)<S_(y)(X′, η_(n))<M⁺(n)

Subsequently, by suitably carrying out transpositions, it can be reduced into an inequality system in X′, similar to inequality system (19). Thus, in a manner similar to the case of Embodiment 1, it is possible to obtain the set of most robust parameter values X′_(s).

In the above-mentioned identifuncation, the signal x may not necessarily be used. Also, the identifuncation may be executed for not only one response of the object system to one signal, but also responses to some different signals. Also, by utilizing the above-mentioned generating convex hull for the suiting, it is possible to provide a more robust suiting region against measurement errors and some noises. Further, in the case of an analog system, it is possible to absorb quantumization errors in the generating convex hull. Yet further, it is possible to apply this embodiment to a design of a system in a time domain by setting a desired time-behavior as the suiting region.

As described above, even when deviations, which have a level corresponding to the robust radius, is introduced into the identified parameter values obtained in accordance with this mbodiment, a behavior of a model still remains within the suiting region. The identified result is most robust against the parameter deviations, among models representing the object system. Thus, in a case where a modeling of an object system rather than an identifuncation by a fixed-point digital filter is merely carried out without taking a degree of the system into consideration, it is possible to obtain a modeling, which is free from deterioration resulting from the quantumization of coeffuncients and input signals, under a condition that a robust radius is greater than a level of quantumization noises, similar to the above-mentioned embodiment.

The obtained robust radius represents a lower limit of the spread of the identifuncation results. Namely, the spread of the identifuncation results is at least larger than the robust radius. Accordingly, the robust radius can be utilized for estimating the precision of the identifuncation. Also, it is possible to obtain an upper limit of the spread of the identifuncation results by forming a convex hull to include the solution region, for example, by calculating the maximums and the minimums of the parameter values.

INDUSTRIAL UTILITY

As is apparent from the foregoing, according to a suiting method and a suiting apparatus of the present invention, on various utilizable problems in industry, such as construction or form of physical systems, identifuncation of systems, display of graphs, data interpolation, signal prediction modeling, pattern recognition, and so on, it is possible to handle input data with a certain spread, which was conventionally diffuncult. Also, by selecting more or most robust solutions from the parameter values of the suiting function, it is possible to construct a higher-robust system against the deviations. 

What is claimed is:
 1. A computer system that determines a set of parameter values of a physical system exhibiting physical characteristics, comprising: an input device that inputs physical characteristics; a memory device that stores a suiting region corresponding to an allowable range of said physical characteristics; a setting device that sets said suiting region, input via said input device, in said memory device; and a processor/controller that executes a calculation to find a suiting function, contained in said suiting region, by solving simultaneous inequalities, said processor/controller including a first setting system that sets said suitng function as a family of functions with parameters corresponding to the parameters of said physical system, a second setting system that sets a distance from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region, and a solution system that solves said siutaneous inequalities so that said distance is maximized thereby determining a robust set of parameter values of said physical system.
 2. A computer system as set forth in claim 1, wherein said processor/controller further comprises a finder that finds approximate distances from the boundary of the solution region of said parameters to parameter values of said physical system.
 3. A computer system as set forth in claim 1, wherein said suiting function comprises a transfer function of an engineering system.
 4. A computer system as set forth in claim 1, wherein said suiting function comprises a recursion relation of an engineering system.
 5. A computer system as set forth in claim 1, further comprising an output device including: a suiting-region-display that displays said suiting region; and a suiting-function-display that displays said suiting function over said suiting region displayed by said suiting-region-display.
 6. A computer system that selects a set of parameter values of a physical system exhibiting physical characteristics, comprising: an input device that inputs a family of suiting functions with parameters corresponding to the parameters of said physical system, said family of suiting functions being contained in a suiting region corresponding to an allowable range of physical characteristics; a memory device that stores said family of suiting functions; and a process/controller that selects a suiting function from said family of suiting functions by solving simultaneous inequalities, said processor/controller further comprising a distance setting system that sets a distance from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region, a solution system that solves said simultaneous inequalities so that said distance is maximized, thereby determining robust a robust set of parameter values of said physical system.
 7. A suiting process performed by a computer system, said comnputer system comprising an input device, a memory device and a processor/controller, to determine a set of parameter values of a physical system exhitibing physical characteristcs, said process comprising: inputting physical characteristics via said input device; setting a suiting region, corresponding to an allowable range of said input physical characteristics, in said memory device; and finding a suiting function, contained in said suiting region, by solving simultaneous inequalities using said processor/controller by: setting said suiting function as a family of functions with parameters corresponding to the parameter of said physical system; setting a distance from a bounds of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region; and solving said simultaneous inequalities so that said distance is maximized, hereby determining a robust set of parmeter values of said physical system.
 8. A suiting process as set forth in claim 7, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and a fourth engineering system represented by a differential equation in time domain.
 9. A suiting process performed in a computer system that includes an input device, a memory device, an output display device and a processor/contoller, to determine and display a set of parameter values of a physical system exhibiting physical characteristics, said process comprising: inputting physical characteristics via said input device; setting a suiting region, corresponding to an allowable range of said physical characteristics, in said memory device; displaying said suiting region on said output display device; finding a suiting function, contained in said suiting region, by solving simultaneous inequalities using said processor/controller by: setting said suiting function as a family of functions with parameters corresponding the parameters of said physical system; settimg a distance from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region; and solving said simultaneous inequalities so that said distance is maximized, thereby determining a robust set of parameter values of said physical system; and displaying said suiting function on said suiting region displayed on said output display device.
 10. A suiting process as set forth in claim 9, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and a fourth engineering system represented by a differential equation in time domain.
 11. A suiting processing performed in a computer system including an input device, a memory device, an output display device and a processor/controller that select a set of parameter values of a physical system exhibiting physical characteristics, said process comprising: inputting a family of suiting functions with parameters corresponding to the parameters of said physical system, via said input device, said family of suiting functions being contained in a suiting region corresponding to an allowable range of physical characteristics; storing said family of suiting functions in said memory device; and selecting a suiting function from said family of suiting functions by solving simultaneous inequalities, using said processor/controller, said selecting further comprising: setting a distance from a boundary of a solution region of the parametes of said suiting function, corresponding to a boundary of said suiting region; and solving said simultaneous inequalites so that said distance is maximized, thereby selecting a robust set of parameter values of said physical system.
 12. A suiting process as set forth in claim 11, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and engineering system represented by a differential equation in time domain.
 13. A memory medium that stores a suiting program executed in a computer system including an input device, a memory device, and a processor/controller to determine a set of parameter values of a physical system exhibiting physical characterstics, said program comprising: inputting physical characteristics via said input device; setting a suiting region, corresponding to an allowable range of said input physical characteristics, in said memory device; and finding a suiting function, contained in said suiting region, by solving simultaneous inequalities, using said processor/controller by: setting said suiting function as a family function with parameters corresponding to the parameter values of said physical system; setting a distance from a boundary of a solution region of the parameter of said function, corresponding to a boundary of said suiting region; and solving said simultaneous inequalities so that said distance is maximized, thereby determining a robust set of parameter values of said physical system.
 14. A memory medium as set forth in claim 13, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and a fourth engineering system represented by a differential equation in time domain.
 15. A memory mediun that stores a suiting program executed in a computer system including an at device, a memory device, an output display device and a processor/controller, to determine and display a set of parameter values of a physical system exhibiting physical characteistics, said program comprising: inputting physical characteristics via said input device; seting a suiting region, corresponding to an allowable range of said input physical characteristics, in said memory device, displaying said suiting region on said output display device; finding a suiting function, contined in said suiting region, by solving simultaneous inequalities, using said processor/controller by: setting said suiting function as a family of functions with parameters corresponding to the parameters values of said physical system; setting a distance from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region; and solving said simultaneous inequalities so that said distance is maximized, thereby determining robust a robust set of parameter values of said physical system; and displaying said suiting function on said suiting region displayed on said output display device.
 16. A memory medium as set forth in claim 15, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and a fourth engineering system represented by a differential equation in time domain.
 17. A memory medium that stores a suiting program executed in a computer system including an input device, a memory device and a processor/controller that select a set of parameter values of a physical system exhibiting physical characteristics, said program comprising: inputting a family of suiting functions with parameters corresponding to the parameters of said physical system, via said input device, said family of suiting functions being contained in a suiting region corresponding to an allowable range of physical characteristics; storing said family of suiting functions in said memory device; and selecting a suiting function from said family of suiting functions by solving simultaneous inequalities using said processor/controller, wherein said selecting of said suiting function comprises: setting a distnce from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region; and solving said simultaneous inequalities so that said distnce is maximzed, thereby selecting robust a robust set of parameter vahles of said physical system.
 18. A memory medium as set forth in claim 17, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and a fourth engineering system represented by a differential equation in time domain.
 19. A process that produces a signal-processing system as a physical system exhibiting specific physical characteristics, using a computer system including an input device, a memory device and a processor/controller, said process comprising: inputing physical characteriscs via said input device; setting a suiting region, corresponding to an allowable range of said input physical characteristics, in said memory vice; finding a suiting function, contained in said suiting region, by solving simultaneous inequalities using said processor/controller, said finding a suiting function comprising: setting said suiting function as a family of functions with parametes corresponding to the parameter of said signal-processing system; setting a distance from a boundary of a solution region of the parameters of said suiting region; and solving said simultaneous inequities so that said distance is maximized, thereby assigning a specific robust set of parameter values to said physical system; and setting parameters, corresponding to the parameter of said found suiting function, in said signal-processing system.
 20. A computer system that determines a set of parameter values of a physical system exhibiting physical characteristics, comprising: an input device that inputs physical characteristics; a memory device that stores a suiting region corresponding to an allowable ramge of said physical characteristics; a setting device that sets said suing region, input via said input device, in said memory device; and a processor/controller that executes a calculation to find a suiting function, contained in said suiting region, by solving simultaneous inequalities, said processor/controller including a first setting system that sets said suiting function as a family of functions with pameters corresponding to the parameters of said physical system, a second setting system that sets a distance from a boundary of a solution region of the parameters of said suiting function, corroding to a boundary of said suiting region, and a solution system that solves said simultaneous inequalities so that said distance becomes larger than a previously-desired level, thereby termining a robust set of parameter values of said physical system.
 21. A computer system as set forth in claim 20, wherein said processor/controller further comprises a finder that finds an approximate distance from the boundary of the solution region of said parameters to parameter values of said physical system.
 22. A computer system as set forth in claim 20, wherein said suiting function comprises a transfer function of an engineering system.
 23. A computer system as set forth in claim 20, wherein said suiting function comprises a recursion relation of an engineering system.
 24. A computer system as set forth in claim 20, frther comprising an output device including: a suiting-region-display that displays said suiting region; and a suiting-function-display that displays said suiting function over said suiting region displayed by said suiting-region-display.
 25. A computer system that selects a set of parameter values of a physical system exhibiting physical characteristics, comprising: an input device that inputs a family of suting functions with parameters corresponding to the parameters of said physical system, said family of suiting functons being contained in a suiting region corrsponding to an allowable range of physical characteristics; a memory device that stores said family of suiting functions; and a processor/controller that selects a suiting function from said family of suiting functions by solving simultaneous inequalities, said processor/controller inluding a distance setting system that sets a distance from a boundary of a solution region of the pareters of said suiting function, corresponding to a boundary of said suiting region, a solution system that solves said simultaneous inequalities so that said distance becomes larger than a previously-desired level, thereby determining a robust set of parameter values of said physical system.
 26. A suiting process that is performed in a computer system including an input device, a memory device and a processor/controller that detee a set of parameter values of a physical system exlutiting physical characterstics, said process comprising: inputting physical characteristics via said input device; setting a suiting region, corresponding to an allowable range of said input physical characteristics, in said memory device; and finding a suiting function, contained in said suiting region, by solving simultaneous inequalities, using said processor/controller by: setting said suiting function as a family of functions with parameters corresponding to te parameter of said physical system; setting a distance from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region; and solving said simultaneous inequalities so that said distance becomes larger than a previously-desired level, thereby determining a robust set of parameter values of said physical system.
 27. A suiting process as set forth in claim 26, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and a fourth engineering system represented by a differential equation in time domain.
 28. A suiting process that is performed in a computer system including an input device, a memory device, an output display device and a processor/controller tiat determine and display a set of parameter values of a physical system exhibiting physical characteristics, said process comprising: inputting physical characteristics via said input device; setting a suiting region, corresponding to an allowable range of said input physical characteristics, in said memory device; displaying said suiting region on said output display device; finding a suiting function, contained in said suiting region, by solving simultaneous inequalities, using said processor/controller by: setting said suiting function as a family of functions with parameters corresponding the parameters of said physical system; setting a distance from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region; and solving said simultaneous inequalities so that said distance becomes larger than a previously desired level, thereby determining a set of robust parameter values of said physical system; and displaying said suiting function on said suiting region displayed on said output display device.
 29. A suiting process as set forth in claim 28, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and a fourth engineering system represented by a differential equation in time domain.
 30. A suiting process that is performed in a computer system including an input device, a memory device, an output display device and a processor/controller that select a set of parameter values of a physical system exhibiting physical characteristics, said process comprising: inputting a family of suiting functions with parameters corresponding to the parameters of said physical system, via said input device, said family of suiting functions being contained in a suiting region corresponding to an allowable range of physical characteristics; storing said family of suiting functions in said memory device; and selecting a suiting function from said family of suiting functions by solving simultaneous inequalities using said processor/controller by: setting a distance from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region; and solving said simultaneous inequalites so that said distance becomes larger than a previously-desired level thereby selecting a robust set of parameter values of said physical system.
 31. A suiting process as set forth in claim 30, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and a fourth engineering system represented by a differential equation in time domain.
 32. A memory medium that stores a suiting program executed in a computer system including an input device, a memory device and a processor/controller that determine a set of parameter values of a physical system exhibiting physical characteristics, said program comprising: inputting physical characteristics via said input device; setting a suiting region, corresponding to an allowable range of said input physical characteristics, in said memory device; and finding a suiting function, contained in said suiting region, by solving simultaneous inequalities using said processor/controller by: setting said suiting function as a family of functions with parameters corresponding to the parameter of said physical system; setting a distance from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region; and solving said simultaneous inequalities so that said distance becomes larger than a previously-desired level, thereby determining a robust set of paramet values of said physical system.
 33. A memory medium as set forth in claim 32, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and a fourth engineering system represented by a differential equation in time domain.
 34. A memory medium that stores a suiting program executed in a computer system including an iput device, a memory device, an output display device and a processor/controller that determine and display a set of parameter values of a physical system exhibiting physical chateristics, said program comprising: inputting physical characterstics via said input device; setting a suiting region, corresponding to an allowable range of said input physical characteristics, in said memory device; displaying said suiting region on said ontput display device; finding a suiting function, contained in said suiting region, by solving simultaneous inequalities using said processor/controller by: setting said suiting function as a family of functions with parameters corresponding the parameters of said physical system; setting a disnce from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region; and solving said simultanmeous inequalities so that said distance becomes larger than a previous-desired level, thereby determining a robust set of parameter values of said physical system; and displaying said suiting function on said suiting region displayed on said output display device.
 35. A memory medium as set forth in claim 34, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and a fourth engineering system represented by a differential equation in time domain.
 36. A memory medium that stores a suiting program, executed in a computer system including an input device, a memory device and a processor/controller that selects a set of parameter values of a physical system exhibiting physical charteristics, said program comprising: inputting a family of suiting functions with parameters corresponding to the parameters of said physical system via said input device, said family of suiting functions being contained in a suiting region corresponding to an allowable range of physical characteristics; storing and family of suiting functions in said memory device; and selecting a suiting function from said family of suiting functions by solving simultaneous inequaities using said processor/controller by: setting a distance from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region; and solving said simultaneous inequalties so that said distance becomes larger than a previously-desired level, thereby selecting a robust set of parameter values of said physical system.
 37. A memory medium as set forth in claim 36, wherein said physical system exhibiting the physical characteristics comprises one of a first engineering system represented by a transfer function of an analog system in frequency domain, a second engineering system represented by a transfer function of a digital system in frequency domain, a third engineering system represented by a recursion relation in time domain, and a fourth engineering system represented by a differential equation in time domain.
 38. A process that produces a signal-processing system as a physical system exhibiting specific physical characeristics, using a computer system including an input device, a memory device and a processor/controller, said process comprising: inputting physical characteristics via said input device; setting a suiting region, corresonding to an allowable range of said input physical characteristics, in said memory device; finding a suiting function, contained in said suiting region, by solving simultaneous inequalities, using said process/controller; and setting parameters, corresponding to the parameter of said found suiting function, in said signal-processing system, said finding a suiting region further comprises: setting said suiting function as a family of functions with parameters corresponding to the parameter of said signal-processing system; setting a distance from a boundary of a solution region of the parameters of said suting function, corresponding to a boundary of said suting region; and solving said simultaneous inequalities so that said distance becomes larger than a previous-desired level, threrby assigning a specific robust set of parameter values to said physical system.
 39. A computer implemented tool that determines a set of parameter values of a physical system exhibiting defined physical characteristics, comprising: an input device through which the defined physical characteristics of said physical system are input; a setting device that sets a suiting region, input via said input device; a memory device that stores said suiting region corresponding to an allowable range of said desired physical characteristics; and a processor/controller that calculates a suiting function contained in said suiting region by solving simultaneous inequalities, said processor/controller further comprising; a first setting system that sets said suiting function as a family of functions with parameters corresponding to the parameters of said physical system, a second setting system that sets a distance from a boundary of a solution region of the parameters of said suiting function, corresponding to a boundary of said suiting region, and a solution system that solves said simultaneous inequalities so that said distance is maximized, thereby determining a robust set of parameter values of said physical system.
 40. A determining device to determine a suiting function possessing a robust set of parameter values corresponding to a set of desired physical characteristics, said determining device comprising: a setting device that sets a suiting region containing one or more sets of parameter values, said suiting region being input through an input device connected with a computer system; a memory that stores said suiting region corresponding to an allowable reage of the desired physical characteristics, said desired physical characteristics bing input through said input device connected with said computer system; and a processor/controller that determines at least one suiting function possessing a robust set of parameter values contained in said suiting region by solving simultaneous inequalities, said processor/controller further comprising; a first setting device that sets said suiting function as a family of functions having parameters, a second setting device that sets a distance measured form a boundary of a solution region of said set of parameters, corresponding to a boundary of said suiting region, and a solution device that solves said simultaneous inequalities so that said distance becomes larger than a previously-determined distance. 